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ABSTRACT 


Approximately 900 unknown atomic energy levels were predicted by 

extrapolation along the helium through sodium isoelectronic Sequences. 

The extrapolations, based on well known regularities in atomic spectra, 
extend beyond the range of known values providing predictions in highly 
ionized atoms. The predicted energy levels are presented, along with the 
known values, in tabular form. In addition, as an aid to sSpectroscopists, 
116 transitions are Listed with known and predicted wavelengths. Since 
the majority of the energy level predictions are in highly ionized atoms, 
most of the predicted wavelengths fall in the vacuum ultraviolet region 


of the spectrum. 
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I. INTRODUCTION 


Identified atomic energy levels provide an indispensable tool for 
scientists in many fields. Scientists in astrophysics, plasma physics, 
physical chemistry, atomic physics and nuclear physics as well as many 
related fields rely heavily on the proper identification of atomic energy 
levels and the corresponding transition values. 

The first decisive step towards a systematic description of atomic 
Spectra was taken by Rydberg, who published his findings in 1890. In the 
80 years following, great progress has been made in the identification 
and classification of atomic energy levels. An excellent selection of papers 
of historical interest may be found in a recent book by Hindmarsh [1]. 

The first major attempt to compile the identified levels resulted in the 
1932 publication of Bacher and Goudsmit [2]. Work continued in spectral 
analysis and soon the number of identified levels increased many fold. A 
more recent compilation of levels was completed by Moore [487] in 1958.° 

Although tens of thousands of levels are now identified, great gaps 
exist in our knowledge. Some spectra are incompletely analyzed and others 
have no identified levels at all. A logical extension of our knowledge 
would be to fill these gaps in order to provide this information for all 
users. 

The indispensable means of energy level identification is analysis of 
experimentally obtained atomic spectra. The experiments on highly 
ionized atoms present difficulties which, in some cases, are almost 
impossible to overcome. In any event, it is helpful to the spectroscopist 


to have an estimate of the transition wavelength to plan optimum use of 


Use of a previously established reference system precludes numbering 
of references in order of occurrence. 











the experimental equipment and to aid in line identification. This esti- 
mate or prediction may be obtained by two methods. One approach would be 
to calculate the energy levels using, for example, the Hartree-Fock method. 
A thorough treatment of this method is given by Slater [3]. The Hartree- 
Fock method is extremely complex for all but the most simple atoms and 
yields results of insufficient accuracy. If exact values could be obtained 
by this method, the complexity of calculations would not be a deterrent. 

The second means of energy level prediction is that of extrapolation or 
interpolation among known values. The validity of this procedure is based 
on well known regularities in atomic Spectra. Extrapolations are relatively 
easy to do in routine fashion but have the disadvantage that extensions can 
be made only along levels which have previously identified values. The 
remainder of this thesis is devoted to the details of this method and results 


obtained. 


II. THEORY 


This section is not indended to be a comprehensive discussion of atomic 
physics. The purpose is to review those concepts which apply directly to the 
extrapolation procedures used. A more detailed discussion of atomic structure 
can be found in many texts, including Herzberg [4], Kuhn [5], White [6] 
and Edlén [7]. 

A. ISOELECTRONIC SEQUENCES 

The name isoelectronic Sequence is used to describe a sequence of an 
atom and ions having the same number of electrons. The sequence may 
Start at any atom in the periodic table; for example, lithium. The neutral 
lithium atom with three electrons is designated Li I and its spectrum is 
frequently called the arc spectrum. The next heavier atom in the periodic 


table is beryllium. Since beryllium has four electrons, the removal of 








one electron will leave the sam number as Li I. This beryllium ion is 
designated Be II and its spectrum is called the first spark spectrum. The 
next heavier atom is boron. To make boron iSoelectronic with Li I and Be II 
two electrons must be removed. This is designated B III and its spectrum 

is called the second spark spectrum. The isoelectronic sequence continues 
in the same manner throughout the periodic table. Since the sequence 
started with the neutral lithium atom it is called the lithium isoelectronic 
sequence. The roman numeral following the element is called the spectrum 


number. 


B. SPECTRAL NOTATION 
lL. One-Electron Systems 

Quantum theory provides a means for specifying the external state 
of the atom. The quantum numbers n (principal), £ (orbital angular 
momentum), s (spin) and j (total angular momentum) are used to specify the 
possible states. The quantity (2s + 1) is called the multiplicity and 
represents the number of values of j. It is common to use letters to 
specify the value of £. These letters correspond to the numerical values: 
s=0, p=1,d=2, £ = 3, g = 4 and down the alphabet skipping the letter 
j. The quantum numbers are not independent but are related as follows: 

nea= 1, 2, 3, 


iis= Ome. 2... .0Ge- 1) 


jzhLrs =h+% but must be > 0 


A typical notation is: 


Which means that: 








mu Lempin. clive 2 
L=1 
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2) Muiltie lectrome systems 


In the multielectron system each electron is described by a set of 
quantum numbers as discussed above. It is necessary to combine or couple 
the angular momenta described by the quantum numbers to form a representation 
of the entire system. Different assumptions about the manner of coupling 
lead to several coupling schemes. Few atoms satisfy completely the 
conditions for any one coupling scheme, so the one chosen represents only 
an approximation of the true coupling. The coupling schemes are labeled in 
a way which denotes the least important interaction among the several angular 
momenta. 

a. L- 5S Coupling 

L - S (or Russell-Saunders) coupling is used almost exclusively 
in this paper and, therefore, will be stressed, In the vector model of. 
the atom, L - S coupling describes the case where the individual spins 
couple strongly among themselves, as do the orbital angular momenta. The 
spin-orbit interaction of the resultants is much weaker. Such coupling 
arises from the predominance of electrostatic over magnetic interactions, 
which is the case in light atoms. For this situation, the same notation 
can be used as for one-electron systems except that the capital letters 5, 
ieeand J refer to the resultant of the individual ee cence The actual 
process of determining the values of S, L, and J for a given electron 
configuration is too lengthly to present here. When equivalent electrons 
are involved, the Pauli exclusion principle limits the number of possible 


states. Let it suffice to give a typical notation example: 


LO}: 








Lee DES 2p* =P 
This represents a system of six electrons with the configuration Se Uasts O ges 
The state of the atom is described by the term = Pa which means that: 

29 + l= 300r%> = 71 

L=l 

J = 0 
This is only one of five possible terms Gps. 2 eed a tee *$,) arising 
Beom this configuration. 

b. Other Coupling 
If the spin-orbit interactions are large compared with the 
electrostatic interactions, as in heavy atoms, j-j coupling becomes 
dominant. Numerous schemes fall in the region intermediate between L - S 
and j-j coupling. Among these are those labeled j - K and L - K, where 
K is the quantum number of the atom's angular momentum (exclusive of the 
spin of the outermost electron). 
3. Parent Terms 
In most configurations, it is impossible to describe the state of 

a multielectron system uniquely with a single term. In such cases it is 
necessary to designate the parent term. As an example, consider the 
nitrogen atom, with seven electrons. If an electron is removed to form 
the ion, the ground configuration of the ion is iee2orep If an electron 
is added now to the parent ion in this configuration, it might go into a 
configuration ls 2s*2p° 3p. Since, however, the ion could have been in a 
state described by See = Pe Pa en or > Diss the state of the atom will 
be different (the energy levels different) depending on the state of the 
parent ion. We describe the state of the atom by including the parent 


term, as ie 2s 2p, Poe eo or lice ce pe@ Ss, ) op “p 
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ne 


Jide 








4. Parity 
If the wave function is reflected through the origin, the function 
is found to be either unchanged or changed in sign only. The function is 
accordingly called even or odd. This property, which is defined for every 
Meemic State, is parity. The parity can be determined from the configuration. 
If the sum of all £'s is even or odd, the parity is correspondingly even 
or odd. Odd parity is designated by a superscript ° on the term, for 


example we 


C. TERM RELATIONSHIPS 

Since members of an isoelectronic Sequence contain the same number of 
electrons, we expect the term systems to be identical except for the 
numerical values of the terms. Thus, it should be possible to predict 
values along a sequence from the known values if a Suitable extrapolation 
formula can be found. 

lL. One-Electron Systems 

One-electron systems are called hydrogenic since the electron 

configuration resembles that of hydrogen. The relativistic quantum treat- 
ment of the one-electron system leads to the following expression for the 


term energy T: 


_ RchZ? , Rcho@ zt ( eae 2) _ Reho” z* 
x 
2 


je a 7 2 3 L + 4n 3 
nN n nN (a) 


(ais +1) - £02 +1) - s(s +1 ; 
2) aia tage) OBS 2 I) 


where: 


ps) 
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Rydberg constant 


velocity of light 


Q 
Wl 


Planck's constant 


2 
Il 


fine structure constant 


R 
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The first term in equation (1) is the Bohr expression for the term energy; 
the second term, the relativity correction; and the third term the spin- 
erbit interaction correction. For each term im the hydrogen isoelectronic 
sequence the only variable in equation (1) is Z. Quantum theory provides 
exactly the term energy values for the hydrogen isoelectronic sequence; 
therefore, the hydrogen sequence is of little interest for extrapolation, 
but equation (1) is quite important in the application to multielectron 
systems. 
2. Multielectron Systems 

In multielectron systems it is necessary to introduce the concept 
of effective nuclear charge, Z - 0, where 0 is the screening constant. 
The screening constant accounts for the fact that in multielectron systems 
the outer electrons are not subjected to full electrostatic attraction by 
the nucleus; part of the nuclear charge is screened by the inner core 
electrons. Substituting Z - 0 for Z in equation (1) gives a fourth degree 
polynomial in Z, assuming O is independent of Z. Actually, Oo is only 
approximately independent of Z but approaches true constancy for large Z. 

3. Energy Level versus Term Value 

The previous discussion applies to term value, which represents 
the energy of the system with respect to the ionization limit. The ioni- 
zation energy is assigned a value of zero and the term value is measured 
down from the ionization limit. The energy level corresponding to the term 
value is customarily measured with respect to the ground state. The ground 
state is assigned a value of zero and the energy level is the value above 
the ground state. Thus the energy level is merely the ionization energy 
minus the term value. Since the ionization energy is expected to demon- 
strate a Z dependence, it should be possible to express the energy level 


by a fourth degree polynomial in Z with appropriate changes in coefficients. 


ils 








[ee RELATIVE TERM VALUES 

Quantum mechanics provides formulae giving the relative positions of 
terms for a given configuration. This knowledge is very useful as a guide 
in empir ical analysis. The following discussion applies only to L-S 
coupling. 

1. Terms of Different L and S$ 

A useful rule is that the lowest level of a given confizuration of 

equivalent electrons (same n and 1) is that with the largest value of S, and, 
if there are several of these, that with the largest value of L. This 
relationship, known as Hund‘s rule, is generally confirmed by observations 
on ground configurations, but many exceptions exist in excited configurations. 
The important fact to note is that regularities do exist which can be used 
as an aid in extrapolation. It is not necessary to rely solely on theory 
in extrapolation, however, as known values should provide the proper 
relationships. 

fee Multiplet Structure 

A multiplet is a set of levels charactérized by the same values of 

L and S, but differing in values of J. The effect giving rise to multiplet 
structure is the spin-orbit term in equation (1). A doublet has two 
allowed J values, a triplet has three allowed J values, etc. If the energy 
of the levels increases with increasing J, the multiplet is described as 
normal. This is usually the case if the unfilled sub-shell is less than 
half filled. If the energy of the levels decreases with increasing J, the 
multiplet in inverted. This is uSually the case if the sub-shell is more than 
half filled. For normal multiplets, the multiplet spacing tends to obey 
an interval rule which was first found empirically by Landé. It states that 
in a normal multiplet the differences between adjacent levels are in the 


ratio of their J values, where for each interval the higher of the two 


14. 








J values is to be taken. Thus, the energy differences “Ps - ape and 
“P, - “Py are in the ratio 2:1. The Landé interval rule can be used to 


obtain the energy level of a multiplet member if the other members are known. 


E. PERTURBATIONS 

The regular arrangement of term values is sometimes found to be 
disturbed. Deviations from the position enced by simple theory are 
called perturbations. One type of perturbation which is well defined is 
the so-called configuration interaction. Under certain conditions, terms 
belonging to different electron configurations will perturb each other. 
This occurs when both electron configurations are of the same parity and 
both terms have the same J value. In addition, in L-S coupling, observation 
shows that the greatest effect is to be expected when the two terms have the 


same L and S values. 


F, TRANSITIONS 

The term value or energy level is not an experimentally measureable 
quantity. However, the energy difference between energy levels may be 
obtained experimentally by spectroscopic measurements. The change from one 
level to another is called a transition. Transitions cannot occur arbitrarily 
between levels but are restricted in dipole radiation. These restrictions 


ermeselection rules are as follows: 


Ag = + 1 (Parity Change) 

Aj = 0, +1 but J=0- J = 0 is forbidden 
1 =0, +1 

AS = 0 


The last two conditions are valid only to the extent that L-S coupling is a 
valid approximation. As we move further into intermediate coupling these 
rules are violated with greater frequency. A transition with AS 70 es 


Si 








called an intercombination transition. These are quite weak in the 

spectra of light elements but become fairly strong in the heavy elements. 
The condition for AJ holds in any coupling. The selection rules forbidding 
dipole transitions between terms of the same parity holds for all 


coupling schemes. 


III. EXTRAPOLATION PROCEDURES 


A. INPUT DATA 
1. Sources 
Over the past several years Professor R. L. Kelly has conducted 

an extensive publications search collecting energy level data. The 
starting point was Moore [487] but hundreds of other references were used. 
This compilation continues today, as pertinent articles appear continuously 
in scientific journals. The data are stored in punch card format, with a 
card made for each energy level for which a value is known. Each card 
contains the following data: 

Element 

Spectrum Number 

Atomic Number 

Configuration and Term Description 

Energy Level Value 

Uncertainty Indicator 

Reference Number 

we iscelectronic Sequence Term Listing 
The initial step requires the organization of the energy level data 

in a manner which simplifies the selection of sequences to be investigated. 
A short program was written which lists the energy level cards for each 
term of the isoelectronic sequence. If the cards are input in ascending 


Ge. 








order within the iso@lectronic sequence; the Jistine by termm 1s also in 
ascending order. This listing shows which terms have many identified 
levels and which terms have few. 
B. EXTENT OF INVESTIGATION 

it was obvious at the outset that the time available would Limit the 
isoelectronic sequences that could be examined. It was also necessary to 
limit the terms extrapolated within each isoelectronic sequence. 

1. JIsoelectronic Sequences 

Because the hydrogen isoelectronic sequence is known completely, 

a logical starting point was the helium isoelectronic Sequence. Investigation 
of ten isoelectronic sequences was selected as a reasonable goal. This 
encompassed all isoelectronic Sequences from helium through sodium. 

2. Terms 

The isoelectronic sequence term listing was used to select terms to 

be extrapolated. (It should be noted here that some predictions are 
actually interpolations but the term extrapolation will be used in all 
cases.) Obviously, five known energy level values are required to obtain 
a unique fourth degree fit. In most cases, then, only those terms with 
five or more known energy levels were extrapolated. (In rare cases 
extrapolations were performed where three or four energy levels were known. 
In these cases only a second degree fit was used.) A further reduction 
was to consider only terms with configurations having all electrons with 


ns 5, 


C. EXTRAPOLATION PROGRAM 
The program used to perform the extrapolations was a standard least- 
squares polynomial fit routine coupled with a plot routine. Input to the 


program consists of the atomic number and energy level value of the known 


War 








isoelectronic sequence members. The least-squares fit routine computes 
the polynomial coefficients and all energy level values up to Z = 30, the 
upper limit selected. The plot routine plots the known input points and 
Superimposes over these the polynomial which gave the best fit. The 
program was restricted to the first through fourth degrees. 
D. EXTRAPOLATION DECISIONS 
ime leeree Used 

The first thing examined was the sign and magnitude of the coeffi- 
cients. If the coefficient of Z was negative, the fourth degree extra- 
polation was discarded, as theory does not support such a relationship. 
Next, the coefficient of Z° was examined. In rare cases this was also 
negative. In every case in which both second and fourth degree coefficients 
were negative, the polynomial plot indicated a possible first degree 
dependence. This is theoretically possible as the subtraction of the term 
value equation from the ionization energy equation could cause the higher 
order coefficients to vanish. In these cases the first degree fit was used 
for extrapolation. If the coefficient of Z* was negative jebue. che 
coefficient of 2° positive, the second degree fit was used. (Perturbations 
and/or incorrect input values were assumed to be responsible as the fourth 
degree fit is much more sensitive.) If the coefficient of Z* was positive 
it was checked for magnitude. In equation (1), the fourth degree terms 
contain on. a very small number. Thus, the coefficient of Z* should be 
much smaller than the coefficient of Z*. If this were not true, the 
second degree fit was used. Again, ee eee and/or incorrect input 
values could have caused the unduly large Z* coefficient. 

2. Validity of Values 
After selecting the degree the output was examined for curve fit. 


Naturally, the better the fit, the more confidence placed on the extrapolations. 
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Curve fits ranged from good to marginal. If the input values had eight 
Significant figures, good curve fits agreed with input values to about six 
significant figures while marginal fits agreed to about four. The extra- 
polated values were only carried to a number of significant figures which 
agreed with the input values. Relationships other than extrapolations 
were used to verify values In multiplets, the Lande interval rule, when 
applicable, was used as a check. In singlets, comparison was made with 
other singlets or with the corresponding multiplet, if present, to insure 
continuation of relationships established in the input values- Violations 
were treated on an individual basis; in some cases the extrapolation was 
used for prediction, in others, the interval estimation was used. 
3. Extent of Extrapolation 

The number of values obtained from a given extrapolation was some- 
what arbitrary, although the limit at Z = 30 was followed throughout. 
Presuming a constant Derceninee of correct input values, the greater the 
input the more valid the extrapolation. Thus, a general guide followed was 
to predict one value for every three input values. This guide was not 
followed when there were only a few known values or when the fit was 


exceptionally good. 
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IV. RESULTS 


The results are expressed in tables of two forms. Tables I, III, V, 

— . XIX list the terms within each isoelectron sequence for which 
Sserapolations were carried out. The columns present, from) leftstom1ene: 
Element, spectrum number, atomic number, energy level value in reciprocal 
centimeters, and reference. Extrapolated ccs are enclosed in parentheses 
and can be further identified by the reference numbers 374 or 375. Those 
known values which have an uncertainty caused by doubtful identification 
or by lack of connection with the ground state are indicated by asterisks. 
In those cases where the extrapolated values differ significantly from the 
known values, and yet the extrapolation appears valid, the extrapolation 
value was listed along with the known value. Also, in cases where the 
multiplet could be resolved by the extrapolation, the extrapolated value 
was listed along with the unresolved value. All levels are designated by 
L-S notation. Since these tables are computer output, certain deviations 
from standard notation were necessary. Subscripts, superscripts and 

roman numerals were not used; therefore, all numbering falls on the line 
and arabic numbers are used throughout. Also fractions were not used, and 
fractional J values are represented by the next higher whole number. 
Finally, an asterisk is used to represent odd parity. 

Tables II, IV, ..... XX list selected transitions within each isoelectronic 
sequence. The tables are presented as an aid to the spectroscopist who 
deals with transitions rather than energy levels. The ee ae are listed 
down the left margin and the sequence members across the top. Those 
transitions involving at least one predicted level are indicated by paren- 
Eaeses. It is impractical to list all conceivable transitions so represen- 


tative transitions which follow normal Selection rules are listed. 
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V. CONCLUSIONS 


The test of the extrapolation procedures is the determination of how 
well the extrapolated levels compare with those deduced from observed 
transitions. Two recent papers list Some energy levels or transitions 
involving energy levels which had been included in this project. 

Tondello [776] observed transitions in the Si XI and Si XII spectrum 
using a laser produced plasma. To aid in line identification, he also used 
extrapolation techniques to predict transitions. Table XXI compares his 
work with predictions of this project. There is reasonably good agreement 
between the two sets of values, the wavelengths agreeing to within 0.1 
Angstrom. 

Gruzdev [GR69] calculated energy levels of the 2p° and 2p 3s 
configurations in the spectra of the oxygen isoelectronic sequence from 
O I through Fe XIX, using the method of intermediate coupling in the single 
configuration approximation. Table XXII compares his work with extrapolations 
Biethis project. Within the 2p configuration the agreement is poor. The 
agreement within the 2p 3s configuration is much better, with less than 
one percent difference in all cases. Comparison by Gruzdev of his 
calculations against experimental values exhibited the same trend, that is, 
much better agreement in the 2p 3s Configuration . 

Of the 48 values compared against other sources 32 were within one 
percent difference. We conclude that it is possible, using extrapolation 
procedures, to predict energy levels with an accuracy that othe their 
use to first approximation. Certainly a first approximation is better than 


no approximation. 
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